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Abstract. An evolutionary scenario to explain the tran- 
sient nature and short total duration of the X-ray burst of 
SAX J1808.4 - 3658 is proposed. An optical companion of 
the neutron star(a"turn-off" Main - Sequence star) fills its 
Roche lobe at the orbital period {Pqj-]j) ~ 19 hours. Dur- 
ing the initial high mass- transfer phase when the neutron 
star is a persistent X-ray source, the neutron star is spun 
up to a millisecond period. Due to its chemical composi- 
tion gradient, the secondary does not become fully con- 
vective when its mass decreases below 0.3 Mq, hence a 
magnetic braking remains an effective mechanism to re- 
move orbital angular momentum and the system evolves 
with Roche - lobe overflow towards a short orbital period. 
Near an orbital period of two hours the mass transfer rate 
becomes so small (~ IO~^^Af0/yr) that the system can 
not continue to be observed as a persistent X-ray source. 
During further Roche - lobe filling evolution deep mixing 
allows the surface of secondary to become more and more 
helium rich . Since the accreted matter is helium rich, it is 
easy to explain observed short total duration of the burst . 
This evolutionary picture suggest that radio emission can 
be observed only at shorter wavelength's . Our model pre- 
dicts a faster orbital period decay than expected if the 
orbital evolution is driven only by gravitational wave ra- 
diation. 

Key words: binaries: close - binaries: X-ray bursters, X- 
ray pulsars (SAX J1808.4 - 3658) 



1. Introduction 

The transient X-ray burster SAX J1808.4 - 3658 was dis- 
covered in September 1996 with the BeppoSAX Wide 
Field Cameras (in't Zand et al., 1998). Two very bright 
type I X-ray bursts with short decay times ~ 8 s (2 - 8 
keV) separated by 14.61 hours were observed. The double- 
peaked time history of both bursts at high energies sug- 
gests a peak luminosity close to the Eddington limit which 
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implies a distance to this object of 4 kpc (in't Zand et 
al. 1998). The observed peak flux of the steady emission 
was about 2xl0~^erg s~^ cm~^ at 2 - 10 keV. For a 
power-law spectrum and a distance, d = 4 kpc the lu- 
minosity is ~ 6x10^^ erg s~^ and for a bremsstrahlung 
spectrum the luminosity is ~5xl0'^^erg s~^. From these 
values the mean accretion rate is ^ 3xl0^^g s~^ (assum- 
ing that the neutron star radius is equal to 10 km and 
the mass is 1.4 Mq). Recently a transient source XTE 
J1808 - 365, positionally coincident with SAX J1808.4 - 
3658 to within a few arcmin, was detected with the Rossi 
X-ray Timing Explorer . The 2-10 keV X-ray flux (Fx) 
on 1998 April 11 was 1.5xlO~Vg s^^cm'^^ MarshaU, 
1998). Wijnands & van der Klis (1998) discovered a co- 
herent 2.49 millisecond signal from this source and they 
estimated an upper limit on the magnetic field strength B 
of ~ (2 - 6) X 10* Gauss. Gilfanov et al. (1998) found that 
between April 25 and 29 1998 this X-ray pulsar showed 
an abrupt decline of the X-ray luminosity which they in- 
terpreted as a result of centrifugal inhibitation at the lu- 
minosity level of a few x 10'^'^ erg s~^. From this they es- 
timated an upper limit on the magnetic field strength of 
B< 3.5x10^ d(kpc) {M/lAMQ)^/^{Fx/10-^ergs-^cm-^ 
Y'^Rr. Gauss, which for d = 4 kpc and a canonical 
neutron star mass (1.4 Mq) and radius (R= 10 km) is 1.4 
xlO* Gauss. Gilfanov et al. (1998) have estimated that 
the total energy of the 1998 outburst was ~ 7.8x10**^ 
erg (3 - 150 keV), corresponding to an accreted mass of 
~1O~^^M0. Since the time between two ourburst was ~ 
1.5 years the average accretion rate is ^10~^^MQyr~^. 

Chakrabarty & Morgan (1998) found the source is 
a binary with orbital period of 2.01 hours, and a very 
small mass function f(M0)=3.85xlO~^M0. Assuming a 
random distribution of binary inclinations, the a priori 
probability of observing a system with inclination i or 
smaller is ( 1 - cos i). For an assumed 1.35 Mq pulsar , 
the 95%-confidence upper limit on the secondary mass is 
M2<0.14:Mq and for a 2 Mq pulsar this limit is O.ISMq. 
Chakrabarty & Morgan(1998) discussed a possible evo- 
lutionary scenario of this system, and proposed that it 
is closely related to the "black-widow" millisecond radio 
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pulsar( PSR 1957+20, PSR J2051 - 0827) which arc evap- 
orating their companion through irradiation. 

Filippenko & Leonard (1998) have found that Keck 
spectra (range 554-685 nm) of the suspected optical coun- 
terpart (Roche et al. 1998, Giles et al. 1998) of SAX 
J1808. 4-3658 reveal a possible very weak, double-peaked 
Ha emission line. 

In this short note we propose a model of the evolu- 
tionary history for this system which is able to explain 
the main observed properties. 

2. Evolutionary picture 

Orbital period determination for low-mass X-ray binaries 
(LMXB) have been a rather difficult task. Up to now ^ 
20 LMXB with the orbital periods less than one day are 
known ( van Paradijs,1995). If we compare the distribu- 
tion of the orbital periods of LMXB and the systems with 
white dwarfs as accretor (catalysmic variable - CV) then 
we can see that the orbital period distribution for two 
classes (LMXB and CV) is different. For both classes the 
orbital periods show a clustering between 3-10 hours. 
There are no LMXB in the CV period gap of 2 - 3 hours 
but if the orbital period distribution for CV shows sec- 
ond peak between 1-2 hours there was no LMXB be- 
tween orbital periods from 50 min to 3 hours (White& 
Mason, 1985). So SAX J 1808.4-3658 is first transient X- 
ray source discovered in this orbital period range. 

Based on data for soft X-ray transients and persistent 
low-mass X-ray binaries with known distance and orbital 
period van Paradijs (1996) showed that for these systems 
the distributions in a (-Pgrb' ■^) diagram can be under- 
stood with the dwarf nova instability criterion adapted to 
account for X-ray heating of the accretion disk. Soft X-ray 
outbursts can occur if the mass transfer rate, M, is below 
a critical value M*'"T, (which depends on orbital period). 
For neutron star low-mass X-ray binaries King et al (1996) 
have obtained the following formula for the critical mass 

accretion rate M^''^, 
crit 



^crit - 5.0 X 10- 



'M%\tJ^Y/^M^yr- 



(1) 



where Mns is the neutron star mass in units of solar mass. 

The evolutionary history of an LMXB and CV depends 
on two main time scales: (i) The hydrogen burning time 
scale inuc and (ii) the angular momentum loss time scale 
iam by magnetic braking, gravitational radiation and by 
mass loss from the system. 

First case 

If the Roche lobe filling secondary is an uncvolvcd 
main-sequence star then inuc ^ ^am- Near the orbital 
period ~ 3 hours the magnetic braking switches off and 
the secondary loses its Roche lobe contact. Second Roche 
lobe filling takes place near the orbital period of 2 hours 
and mass exchange rate is now « lO""'^" M0/yr. Since this 



mass accretion rate is larger than M""^. the persistent X- 
ray sources with short orbital periods must be observed. 

In the period gap the evolutionary paths for CV and 
LMXB may greatly differ. A neutron star with sufficiently 
weak magnetic fields accreting for a long time from a 
surrounding Keplerian disk can be spun up to millisec- 
ond periods. Van den Heuvel&van Paradijs (1988) pro- 
posed that the evaporation of companion stars (induced 
by MSP) could account the lack of persistent low-mass 
X-ray sources below the period gap. 

Second case 

If the Roche lobe filling secondary is an "turn-off' 
Main-Sequence star thentam"^inuc- For this case accord- 
ing to the results of computations (Tutukov et al 1985, 
Pylyser & Savonije , 1989, Ergma et al , 1998 ) at first 
the mass exchange proceeds on the thermal time scale of 
the secondary, but later M quickly decreases. Due to a 
chemical composition gradient the secondary did not be- 
come fully convective when its mass had decreased be- 
low 0.3 Mq. Magnetic braking is still the most effective 
mechanism to remove orbital angular momentum and the 
system evolves towards a short orbital period via Roche - 
lobe overflow. In the period gap the mass accretion rate 
is very low (^ 10~^°M©/yr) which is less than the criti- 
cal mass accretion rate for system to be persistent X-ray 
source. 

First case allows easily to understand the evolution of 
the LMXB with the orbital period between 3 and 10-12 
hours. 

Since SAX J1808.4 -3658 is transient X-ray source we 
propose that the progenitor SAX J1808.4 - 3658 evolves 
according to second case . Using the evolutionary program 
of Sarna & De Greve(1994,1996) we have calculated an 
evolutionary sequence for a binary with a 1.4Mq neutron 
star and a IMq companion. The secondary fills its Roche 
lobe a -Pqj.^3 = 0.8 days. When the orbital period is two 
hours the M, is 4.7x10^^^ M^yr'^ and the surface hy- 
drogen abundance has been decreased to X = 0.3. As we 
will see later this value is important to explain the short 
duration of the X-ray bursts. 

In Fig.l we have drawn the value of the critical mass 
accretion rate as a function of ^orb' Near an orbital pe- 
riod of two hours the calculated mass accretion rate is 
below the critical value and hence the system can not be 
a persistent X-ray source. 

Chakrabarty & Morgan (1998) proposed that if the X- 
ray bright state lasts long enough (or if the pulsations re- 
main detectable in quiescence) , further X-ray longer-term 
timing observations may yield a detection of orbital pe- 
riod evolution, which would probe effects of gravitational 
radiation, tital interactions and mass loss. If the orbital 
evolution is driven by gravitational wave radiation then 
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Fig. 1. Evolution of the mass accretion rate versus orbital 
period. The critical mass accretion rate ( equation (1)) is also 
shown (dashed line). The arrow shows the position of SAX 
J1808.4 - 3658 



the orbital period will decay according to (Landau & Lif- 
shitz, 1971) 



P, 



1927r, 2tt 



orb 



orb 



P. 



5/3M^Mns Mr 



Mq Mq 'M, 







'1/3(^)5/3(2) 



GMr. 



T= 



- = 4.925 xlO ^s Ms is the mass of secondary 
Ms +Mns is the total mass . If the secondary 



where 

and M 

has a mass of O.O5M0 and Mns=l-4M0 (as discussed by 

Chakrabarty & Morgan, 1998) then P^^^^ 8.1 xlQ'^^. 

In our model, since the orbital evolution is driven by mag- 
netic braking the orbital period decay is faster (^— 4.5 
xlO~i^). So the determination of the value of orbital pe- 
riod decay may give important clue to understanding the 
prehistory of this system. 

2.1. The chemical composition of the secondary 
and the X-ray bursts 

In our scenario the secondary is a helium rich star. There 
are two systems, namely MXB 1916 - 05 (^orb ~ ^^ min) 
and MXB 1820-30 (PQj.b==11.4 min), which according to 
current understanding of the low-mass binary evolution 
with compact object as an accretor, must have a Roche - 
lobe filling helium rich secondary in order to explain its 
short orbital period (Paczynski & Sienkewicz ,1981). 

Computations of X-ray bursts at the surface of accret- 
ing neutron stars have shown (Ayasly&Joss,(AJ) 1982, 
Kudrjashov&Ergma,1983) that a short total duration rj 
of the burst is possible only if the matter is helium rich 
(X<0.2-0.3, where X is the hydrogen abundance by mass) 
before unstable burning starts. It may be possible that 
during the accretion phase the hot CNO hydrogen burn- 
ing takes place in the accreted matter on a timescale ^ 
10^ X {N-p/Nq-^q) s « 10'"' s. If the recurence time be- 
tween two succesive burst is longer that this hot CNO 



hydrogen burning time the hydrogen may be exhausted 
at the bottom of the accreted envelope . Crude estimate 
using the bursts light curve published by in't Zand et al 
(1998) gives rj ^ 14 - 15 s. Let us analyze the observed 
data for the X-ray burst using the AJ numerical models. 
Since the accretion rate (from observations) is a few times 
]^Qi6g g-i lY^Qj^ ^Q may compare it with the AJ model 1 
which has been computed with a similar mass accretion 
rate (Z — 0.01 where Z is the metallicity of the accreted 
matter). In this model the hydrogen abundance, ATu, at 
the base of the freshly accreted matter prior to the flash 
is 0.59. The burst total duration , tj, is 34 s and the 
burst recurence time, irec ~ 4.7 hours . The maximum 
burst luminosity Lmax is 5xlO'^^erg a~^. The AJ model 
3, in which the mass accretion rate is one order magni- 
tude less, has X^ = 0.0004, tree ^ 16.3 hours, t^^ ^ 7 
s and imax = 2.5xl03'^erg s"^ The AJ model 16 with 
the same accretion rate as for model 1, but with Z=0.001, 
has X^ — 0.65, tree ^11.2 hours, rj ^ 38 s and Lmax 
= 7.5xl0'^^erg s"^. From these results it can clearly be 
seen that in order to have a short total duration of the 
burst , the hydrogen must either be exhausted or the ac- 
creted matter must be helium rich in the envelope. In SAX 
J1808.4 - 3658 the time between two observed successive 
bursts is ^ 14 hours, which is much longer than for model 
1 with a similar mass accretion rate. One possible expla- 
nation is that the amount of CNO elements is reduced due 
to spallation, as it was discussed by Bildsten et al. (1992), 
and hence the recurrence time will increase ( compare AJ 
models 1 and 16). It is also possible to speculate that the 
burst recurence time may increase with increase of the ra- 
dius of the neutron star (referee comment). For example 
in AJ model 14 (_Rns=13 km) the recurrence time is ~ 
16 hours. In our analyze we used A J standart model re- 
sults ( i?ns ~ 6.57 km ) since according to in't Zand et al 
(1998) a burst emitting sphere radius is of 8 km ( at the 
distance of 4 kpc) . Although the neutron star radius was 
determined by in't Zand et al using a simple black body 
fit to the burst spectrum and can therefore underestimate 
the radius of the emitting sphere (Titarchuk, 1994). 

So to explain the short total duration of the burst and 
the Eddington luminosity we propose that the accreted 
matter must be helium rich . In order to explain long 
recurrence time one would need either reduced Z value or 
large i?ns- 



2.2. The possibility of detection of the millisecond 
radiopulsar. 

Since the millisecond X-ray pulsar SAX J1808. 4-3658 has 
orbital parameters very close to those of the eclipsing bi- 
nary pulsar system PSR J2051-0827, it has been proposed 
(Chakrabarty & Morgan, 1998, Wijnands & van der Klis, 
1998) that this system may emerge as a radio pulsar after 
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the X-ray state ends. Lets us estimate the optical depth 
for the free - free process (Illarionov&Sunyaev,1975) 



^ff = 



100(M/MEdd)\ A T 



-i.5/f_gTh\- 



lyr 



(3) 



where Afgj j is the Eddington mass accretion rate. Ac- 
cepting M- lO-"M0yr~\ T=6000 (less than the hydro- 
gen ionization temperature), Mns=1.4 Mq,Ms = O.IM0, 
A= 3 cm and ^orb~ ^ hours we can find that tq is ^ 
1, which means that at A < 3 cm it may be possible to 
observe radio emission. For longer wavelenghts (20 or 70 
cm) tq is greater that unity. 

3. Conclusion 

We propose that the progenitor of SAX J1808.4 - 3658 
was a low mass binary with a neutron star and a secondary 
that filled its Roche- lobe at the turn-ofi^ Main-Sequence 
.During the bright X-ray phase the neutron star will spin- 
up to a millisecond period. Due to a chemical composition 
gradient, the secondary does not become fully convective 
when its mass has decreased below 0.3 Mq and magnetic 
braking is therefore still the most effective mechanism to 
remove orbital angular momentum. Near an orbital pe- 
riod of two hours the mass transfer rate is so small {^ 
lO~"'^"'^M0yr~-'^) that the system can not be a persistent 
X-ray source. In our evolutionary scenario, due to deep 
mixing, the secondary is a helium rich star. Since the ac- 
creted matter is helium rich it is easy to expain the ob- 
served short total duration of the burst . It is estimated 
that in this evolutionary picture it may be possible to ob- 
serve radio emission from the pulsar only at short wave 
lenghts, A< 3cm. Our model predicts a faster orbital pe- 
riod decay than in the case where the orbital evolution is 
driven only by the gravitational radiation. 
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